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I .  INTRODUCTION 


Let  y.21  •••!  Independent,  Identically-distributed  random 

variables.  Then  <  Y2  <  <  Y^^,  where  the  Y^'s  are  the  X^^'s  rearranged 

in  order  of  increasing  magnitudes,  are  defined  to  be  the  order  statistics 
corresponding  to  the  original  random  sample.  Order  statistics  find 
immediate  application  in  the  design  and  evaluation  of  logical  structures 
which  make  decisions  based  on  the  relative  values  assume.:  hy  a  set  of  n 
random  variables.  One  particularly  interesting  appil.->  ion  Involves  the 
determination  of  time  windows  for  firing  Impulses  in  a  fuzing  system. 

This  is  the  problem  which  motivated  the  work  on  which  we  are  reporting. 
However,  the  work  is  general  In  nature  and  may  prove  useful  for  other 
applications  in  the  armament  research  and  development  community, 


II.  STATEMENT  OF  THE  PROBLEM 

In  the  particular  problem  which  we  addressed,  a  fuze  contains  N  detonators, 

K  of  which  must  function  within  a  specific  time  span.  Furthermore,  the  second 
detonator  (which  functions  at  time  Y2)  partitions  the  time  span  into  two  sub¬ 
intervals.  The  first  subinterval  [Y2-  Yj]  is  examined  to  determine  if  the 

first  detonator  functioned  within  that  time  segment,  and  the  second  subinterval 
[Y2,  Y2  *  <$2]  is  monitored  to  count  the  number  of  additional  detonators  activated 

during  that  period  of  time.  If,  within  the  time  Interval  [^2  -  Y2  +  62I' 

K  detonators  have  functioned,  then  the  command  to  fire  will  be  initiated; 
otherwise,  it  will  not.  The  times  to  function  for  the  detonators  are  random 
variables  and,  as  such,  can  be  characterized  by  a  cumulative  distribution 
function  P.  Assuming  that  the  time  to  function  of  each  detonator  is  identi¬ 
cally  distributed,  then  the  problem  consists  of  expressing  the  probability 
of  fuzing  as  a  function  of  K,  N,  5^,  62,  and  P. 


III.  SOLUTION 


Let  be  the  time  to  function  of  detonator  i  in  its  operating  environ¬ 
ment.  Then  X^,  X2,  ...,  X^^  are  independent,  Idei  tically-dlstrlbuted  random 
variables;  and  we  can  define  Yj^,  Y2,  ...,  Y^^  to  je  the  order  statistics 
corresponding  to  the  X^'s.  Por  our  problem,  if  Y2  "  ^  “’f®  interested 

in  the  probability  that  Yjj  -  Yj  <  62*  however,  if  ^2  -  Yj^  >  6j,  we  need  to 
determine  the  probability  that  Y|^^j  “  ^2  ^  '^2  *^**“*"^^8  K  +  1  <  N.  Therefore, 
we  need  to  evaluate 


Pr  {warhead  fuzing}  ■  Pr  {Yj  -  Yj  <  Pr  {Yj^  "  ^2  ^  *^2  ^  ^2  ‘  ^1  ^  '*1^ 


+  Pr  {Yj  >  Yj  >  fij}  Pr  '  ''^2  ^  ®2  I  ^2  ■  '^1  ^  • 


(1) 
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which  !•  equal  to 


■’l  ■  -(K-j)  liN-K)l  ft''*"’  ■  ’■'’'“’l'''’ 

•  f(u+v+w)  [1  “  F(u+v+w)]^"^  du  dv  dw  .  (8) 

In  an  analogous  manner  we  can  obtain  the  Joint  probability  density 
function  of  the  1st,  2nd,  and  K+lst  order  statistics;  and,  letting 
u  ■  yjj^.2  "  y2*  "  ^2  '  ^1  ^  *  ^1’  obtain  the  necessary 

probability  for  the  second  term  on  the  right-hand  side  of  Equation  3. 

That  probability  is  equal  to 

4>00  4*^  (S 

*  f(u+v+w)  [1  -  FCu+v+w)]*^"'^‘^  du  dv  dw.  C93 


The  probability  of  fuzing  is  then  Just  the  sum  of  Equation  8  and  Equation  9. 

Appendix  B  contains  a  computer  program  which  evaluates  these  inte¬ 
grals.  In  its  current  form  it  assumes  that  the  distribution  of  the 
times  to  function  of  the  detonators  is  normal;  however,  it  can  be  easily 
modified  to  change  this  distribution.  The  program  requires  as  input  N, 

K,  (Sj,  62,  and  a  (standard  deviation  of  the  assumed  distribution). 

IV.  RESULTS 

For  the  problem  we  addressed,  the  value  for  a  was  specified  to  be 

10"®  seconds.  Figure  1  presents  the  results  of  changing  6j  and  §2  ® 

fuzing  system  in  which  K  is  equal  to  N-1.  In  Figure  2  we  considered  a 

fuze  with  eight  detonators;  and,  keeping  0  equal  to  lO'^,  we  varied  K  as 
well  as  and  62’  course,  given  any  four  of  the  five  input  variables 

(o,  K,  N,  and  62),  we  can  obtain  a  specific  probability  of  warhead 

fuzing  by  parametrically  varying  the  remaining  variable. 


V.  APPLICATION  TO  FUZE  DESIGN 

An  important  consideration  is  to  determine  optimal  6^  and  6^  based  on 

a  specified  0,  K,  and  N  and  a  desired  probability  of  warhead  fuzing.  To  do 
this  several  factors  must  be  considered.  To  maintain  the  reliability  of  the 
fuze  it  is  necessary  that  the  interval  [Y2  "  *5^,  Y2  +  62!  sufficiently 
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FIGURE  1.  PROBABILITY  OF  IMRHEAO  FUZIHG  IM  AN  H-1  OUT  OF  N  FUZING  SYSTEM  (o  = 


FIGURE  2.  PROBABILITY  OF  WARHEAD  FUZING  IN  A  K  OUT  OF  8  FUZING 


large  to  allow  for  the  occurrence  of  K  detonations.  However,  safety  must 

also  be  considered;  and  here  it  is  desirable  that  the  interval 

[Y2  “  ^2  ^  small  as  possible  to  preclude  the  occurrence  of 

spurious  detonations.  These  conflicting  goals  could  be  accomplished  through 
a  solution  to  the  following  optimization  problem; 

Given  0,  K,  and  N,  determine  6^,  62  >  0  for  which 
dj  +  62  is  minimum  subject  to  Jj  +  J2  >  ex  (reference 
Equation  8  and  Equation  9]  where  a  is  the  desired 
probability  of  warhead  fuzing. 

Consideration  of  this  problem  is  beyond  the  scope  of  this  report.  Here, 
we  will  content  ourselves  with  the  determination  of  reasonable  initial  values 
of  6^  and  62  for  use  in  the  computer  algorithm.  Such  values  may  then  be 

adjusted  as  required  to  obtain  the  desired  probability  of  warhead  fuzing. 

For  this  purpose  practical  initial  values  for  6^  and  ^2  would  be  B[Y2  -  Y^] 

and  E[Yj,  -  Y2]  respectively,  where  E[.]  is  the  expected  value  of  a  random 

variable. 

For  the  general  case  we  can  define  Wj^  ■  "  ^R*  ‘^o>'®ideration 

of  the  material  in  Appendix  A  we  obtain  the  probability  density  function 

00 

•  [FCX  +  w^g)  .  F(x)]®“^“^  f(x  +  w^g) 

•  [I  -  F(x  +  dx  (10) 

where  Cj^g  -  (r.j)  i  (s-r-i)  1  (N-sTT  * 

To  determine  "  B[Yg  -  Yj^]  we  must  evaluate 

4.00 

"  ^0  '*RS  ^^*'rS^  ‘^''rS* 

Therefore,  we  can  obtain  EtYj  -  Yj]  by  evaluating 

4.00  «fOO 

B[Wj2]  -  C^2  f,  ^12  f(x  4.  Wj2) 

•  [1  -  P(x  +  dx  dWj2;  (12) 
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and  in  an  analogous  manner  we  can  obtain  E[Y|^  -  Y2]  by  evaluating 


E[W2|(]  ■  Cju  W2K 

•  [FCx  ♦  WjK^  “  *  «2K^ 

•  [1  -  F(x  +  '^'^2K  ■ 

The  pth  quasi-range  of  N  sample  values  is  defined  to  be  the  range  of 
|>)  -  2p)  values,  omitting  the  p  largest  and  the  p  smallest.  For  the  special 
case  in  which  K  ■  N-1,  -  Yj^_j  -  Y2  becomes  the  first  quasi-range;  and 

this  distribution  has  been  investigated  and,  to  some  extent,  tabulated. 

For  example,  H.  Leon  Harter*  has  provided  the  expected  values  of  the  first 
quasi-range  »  ElWj  for  samples  from  a  normal  distribution  with 

mean  zero  and  variance  one  and  for  values  of  N  from  four  to  one  hundred. 
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APPENDIX  A 

DERIVATION  OF  THE  JOINT  PROBABILITY  DENSITY  FUNCTION 


As  shown  by  H.  A.  David,  we  can  obtain  the  joint  probability  density 
function  for  Y^,  Yj,  and  Yjj.  Lot  Xj,  X2,  . . . ,  be  a  random  sample  from 

a  population  with  density  f(x)  ■  F'{x);  and  let  Y^,  Y21 

corresponding  order  statistics.  Then  the  probability  density  function  for 
the  rth  order  statistic  may  be  derived  by  considering  the  following  config¬ 
uration; 


That  Is,  Xj^  <  X  for  r-1  of  the  X^^,  x  <  X^  <  x  +  Ax  for  one  X^^,  and 
X^  >  X  +  Ax  for  the  remaining  n-r  of  the  Xj^.  The  number  of  ways  this 
combination  of  events  can  occur  Is 


and  each  such  way  has  probability 

[F(x)]’^'^  [P(x+Ax)  -  F(x)]^  [1  -  F(x+Ax) 


Therefore,  we  have 

Pr  {x  <  Y^  <  X  ♦  Ax}  -  ^ir:-iyT- (nT)  !’ 

•  [FCx)]’^"^  [F(x+Ax)  -  F(x)]  [1  -  FCx+Ax)]"’^‘  +  O(Ax^)  (Al) 

2  2 

where  0(Ax  )  means  terms  of  order  (Ax)  and  Includes  the  probability 

of  realizations  of  x  <  Y.  <  x  Ax  in  which  more  than  one  X.  Is  In 

r  1 

(x,  X  *  Ax) .  Dividing  both  sides  of  Equation  Al  by  Ax  and  then  letting 

Ax  0,  we  obtain 


fy  CX) 
r 


[FCx)]’"'^  f(x)[l-P(x) 


David,  H.  A.,  Order  Statistics,  John  Wiley  and  Sons,  Now  York,  1970,  Pp,  9-10. 


In  a  similar  manner  we  can  derive  the  Joint  probability  density 
function  of  Y^,  and  Yg! 


r-1 

8-r-l 

n-s 

X 

x+Ax  y 

y+Ay 

and  through  an  analogous  argument  we  obtain 

7  9 

•  lF(y)  -  £(y)  [1  -  FCy)]"'®.  (A3) 

Finally,  for  the  joint  probability  density  function  of  Y^,  Y^,  and  Y^ 


r-1 

I 

s-r-1 

L 

t-8-1 

1 

n-t 

X 

x+Ax  y 

y+Ay  i 

{••■Az 

and 

^Yy.Yj,Yt^*'^'*‘^  "  "(r-13  >  (s-r-1)  I  (t-s-l)  1  (n-t)  1 
•  tF(y)-F(x)]*“’'"^  f(y)  tF(z)-F(y)]^"*'^  fCz)  [l-F(^)]"■^  (A4) 

For  the  case  r-1,  5*2,  and  t-k.  we  obtain 

^YpY2,Yk^''*y'*^  “  ■nrsTrin^kyi 

•  [F{z)-F(y)]''"^  f(z)  (l-P(z)]"-‘'.  (AS) 
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APPENDIX  B 
COMPUTER  PROGRAM 

N«  are  presontlng  here  the  computer  program  which  evaluates  the  desired 
probability  function.  As  noted  In  the  program  comments,  we  have  assumed  that 
the  times  to  function  of  the  detonators  are  normally  distributed  with  mean 
2 

zero  and  variance  a  .  However,  with  Just  a  few  changes  to  the  subroutines, 
a  different  distribution  may  be  assumed.  To  do  this,  all  cards  containing 
'7<K)IMAL'*  in  columns  74  through  79  must  be  replaced  by  others  with  the  approp¬ 
riate  probability  density  function  or  cumulative  distribution  function. 


wuuuwuuwuoouiruuwououovuvu  «  o  u  u  u  uouo 


PRUOMM  OtToN  (INPWT»OUTMUT»TAI«eS»INRUT,T*Pt6»OUTPUT) 


THIS  PROaSAM  OeTKRMXNIt  THK  PROBAIILZTY  THAT  K  OUT  OP  N 
OCTONATORS  pZU  PUHCTION  WITHIN  A  OIVfN  TIMCSPAN,  K  MUST 
■C  ORCATIR  THAN  OR  CflUAL  TO  S*  ANO  K  MUST  SC  LESS  THAN  OR 
eoUAC  TO  N«  THE  OISTRlpuTtON  OP  THEIR  PUHCTIOnInO  IS  ASSUHCO 
TO  II  OAUSIIAN  WITH  A  HIAN  lOUAL  TO  e«0  AND  A  ITANOARO 
DEVIATION  lOUAL  TO  SIONA,  THE  SOLUTION  IS  DERIVED  THROUOH 
THE  USE  OP  ORDER  STATIITICB  AND  IS  OBTAINED  BY  CVALUATINO 
A  TRIPLE  INTEORALi 

IP  THE  NUMBER  OP  DETONATORS  IS  OREATER  THAN  TNIRTVt 
THEN  THE  TOLERANCE  LIMITS  OP  THE  INTEORALS  MUST  BE  REDEPINEOi 
THAT  ISi  THE  VARIABLE  <CRROR«  IN  THE  MAIN  ROUTINE  ANO  THE 
variables  'TOLXtTOLVt  IN  SUBROUTINE  'OtST*  SHOULD  BE  ADUUSTEO. 
with  THE  tolerance  limits  CURRENTLY  IN  THE  PROORAH* 

The  RESULTINO  probabilities  are  CORRECT  TO  TWO  decimal  PCACES. 

INPUT  IS  AS  POLLOWS 


N  ...... ...........  number  op  DETONATORS 

K  . . . . NUMBER  OP  DETONATORS  TO  PUNCTION 

SIOMA  .............  STANDARD  DEVIATION  OP  DISTRIBUTION 

OELTI  . . .  TIMESPAN  (PIRIT  TO  SECOND  OETONATORS) 

OElTS  . .  TIMESPAN  (SECOND  TO  LAST  DETONATORS) 


REAL  NPACTiKZPACTiKSPACT.NKPACTtNKlPACT 
DIMENSION  IERRIB) 

COMMON  /COMI/  N.K.SIOMA.OELTItOELTt 
COMMON  /COMI/  PI«W 
COMMON  /COMA/  I NO 

external  OIST 

data  IERR  /3*(-0)«0«»p(««)/ 

CALL  SYSTEMC  OA.IERRI 
CALL  SYSTEMC  (IlSiIERR) 

WRITE  (AilSO) 

READ  INPUT 

initialize  variables 

5  READ  (B.lOO)  N.K.IIOMA.OILTI.OELTI 
IP  (K.LT.3  .or.  K.OT.N)  00  TO  IS 
IP  (N  .OT.  30)  WHITE  lOtSOO)  N 
IP  (OELTI  .OT.  lO.wSlOMA)  DILT1«10.*SIOMA 
IP  (OELTI  .OT.  lO.MIOMAl  0ClT|w10.*SI0MA 
RCSlwO. 

PEStwO. 

Pla3.;AI99IOS30 

NNaN/S*0 
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uoo  WOU  tfW  u  uuu 


If  (N  .Cfl.  80)  NNnlC 

eRROHaIOt**(«NNI 

Aa«10.*IXOM 

eVALUATC  TMIfLK  INTEOIIAC- 

INOal 

PCSlaSaUANK  lA.I.CRRORtRUM.OIST) 

If  (N  .ea«  R)  00  TO  T 

INOal 

RCSIaSOUANK  (AiltCRRORtRUMtOISTl 
OCTCRMINC  fACrORIAUI 

7  CONTINue 
NfACTal 
KlfACTat 
KSfACTal 
NKfACTat 
NRlfACTal 
00  10  laliN 
NfACTaNfACT*! 


If 

(I 

•  IE. 

(K-2) ) 

K2fACTwK2fACTa| 

If 

<1 

.LE. 

(K-3)  ) 

KSfACTaKSfACT*! 

If 

(I 

.UE. 

(N-K) ) 

NKfACTaNRfACT*! 

If 

(X 

.).E.- 

<N«K-1))  NRlfACTaNKlfACT*! 

10  CONTINUE 

COMPUTE  AND  PRITE  PROOAOICXTIES 
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